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Abstract. Finite index subgroups of the modular group are of great 
arithmetic importance. Farey symbols, introduced by Ravi Kulkarni in 
1991, are a tool for working with these groups. Given such a group F, a 
Farey symbol for F is a certain finite sequence of rational numbers (rep- 
resenting vertices of a fundamental domain of F) together with pairing 
information for the edges between the vertices. They are a compact 
way of encoding the information about the group and they provide a 
simple way to do calculations with the group. For example: calculating 
an independent set of generators and decomposing group elements into 
a word in these generators, finding coset representatives, elliptic points, 
and genus of the group, testing if the group is congruence, etc. In this 
expository article, we will discuss Farey Symbols and explicit algorithms 
for working with them. 



1. Introduction 

Modular forms are certain functions defined on the upper half plane dis- 
playing certain symmetries under the Mobius transformation action of a 
finite-index subgroup of PSL<2(Z). The theory of modular forms has been in 
the central stage of number theory for more than one century and continues 
to be one of its most exciting areas. Working with modular forms requires 
knowing information about their underlying groups. There is a vast liter- 
ature about many aspects of finite index subgroups of the modular groups 
and their relations with other fields such as combinatorics, algebraic curves. 
Interested readers are referred to articles like [ASD71], [Bir94], or a recent 
survey article by the second author on these groups and their modular forms 
[Lon07]. 

Some finite index subgroups of the modular group can be described purely 
by congruence relations, and as such are called congruence subgroups of 
the modular group. These groups are relatively easy to work with, as they 
contain a certain normal subgroup T(N), such that the quotient of the group 
by r(iV) is just a subgroup of PSL2(Z/iVZ). Most computational methods 
for working with modular forms work only for congruence subgroups (cf. 
[Ste07]). Noncongruence subgroups are finite-index subgroups of PSL^Z) 
that cannot be described by congruence relations. Even though they make 
up the majority of finite- index subgroups of PSL/2(Z) there are few tools for 
working with them. One tool that works equally well with both congruence 
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and noncongruence groups is the method of Farey symbols introduced by 
Ravi Kulkarni [Kul91]. In this expository paper we will recast how to use 
Farey symbols and some related computational topics. We will discuss some 
explicit algorithms for working with Farey symbols. The first author has 
implemented a collection of such algorithms into a free SAGE package called 
"KFarey". It should be made clear to the readers that "KFarey" is an 
undergoing project and we will continue to improve the current functions, 
implement other existing algorithms such as [Hsu97, Lan02] , and investigate 
new algorithms to make "KFarey" useful to a wide range of audience. The 
second author would like to thank ICAC 2007 organizers for inviting her to 
attend the conference and give a talk on this topic. 

2. Subgroups of PSL 2 (Z) 

Let SL 2 (Z) be the group of 2 x 2 matrices with integer coefficients and 
determinant 1, and let PSL 2 (Z) = SL2(Z)/{7, — I}. Let M be the upper half 
plane H := {z G C : lm(z) > 0}. Then PSL 2 (Z) acts faithfully on H under 
the action 

az + b 
cz + a 

where 7 = ^ ^ € PSL 2 (Z). The objects of our study will be the finite 

index subgroups of PSL 2 (Z). For example, the standard congruence groups, 
T(N), Ti(N), and r (JV). 

If r is a finite index subgroup of PSL 2 (Z) then the action of T partitions 
Q U {oo} into equivalence classes, where q\ ~ g 2 if q\ = 7g 2 for some 7 G T. 
These equivalence classes {q} are called the cusps of T and the width of 
the cusp {q} is [Stab PgL2 ( Z )(g) : Stabr (<?)]• We say that the level of T is the 
least common multiple of the cusp widths of T. 

Recall that a congruence subgroup of PSL 2 (Z) is a subgroup T that con- 
tains a principal level N congruence subgroup T(N) for some N. If TV is the 
smallest such N such that this is true, then F has level N. 

Definition 1. Let T be a finite index subgroup o/PSL 2 (Z). For our pur- 
poses, a fundamental domain of T is a hyperbolic polygon P on HuQujoo} 
such that: 

(1) If z is in the interior of P and 7 G T, then 72 £ P implies 7 = I. 

(2) For every z € H there is 7 € T such that 72; G P. 

Lemma 1. Let V be a finite index subgroup o/PSL 2 (Z) and P a hyperbolic 
polygon. Suppose P is such that: 

(1) If z is in the interior of P and 7 G T, then 72 G P implies 7 = I. 

(2) For each side e of P, there is 7 G T such that 7 maps e to another 
side of P in an orientation reversing manner. 

Then P is a fundamental domain ofT. 
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Proof. Let the images 7P of P under elements 7 of V be called P-tiles. By 
Condition 1 they cannot overlap. Also, given Condition 1, we only need 
to show that H C [^J 7P. Suppose this is not true. Then there is n G V 

such that r/P has an edge e without a P-tile on the other side. Then rf x e 
is an side of P and the 7 of Condition 2 maps rf x e to another side of P. 
Specifically, 7 _1 maps P to a P-tile adjacent to P across the side e. Then 
r/7 _1 P is a P-tile adjacent to r/P across the side e. Contradiction. □ 



3. Farey Symbols 

3.1. Special Polygons. Farey Symbols were introduced by Ravi Kulkarni 
in 1991 [Kul91] as a compact and efficient way to compute with finite index 
subgroups of PSL2(Z). The idea is to describe the group by a fundamental 
domain with vertices at certain rational numbers and certain hyperbolic arcs 
joining these rational numbers. Most of the theory here is summarized from 
[Kul91]. 

If x and y are two points on H U Q then there is a unique circle passing 
through x and y with center on Q. We say the hyperbolic arc joining 
x and y is the arc of this circle contained in EI U Q joining x and y. We 
also say the hyperbolic arc joining x G EI to 00 is the vertical line segment 
{x + ti : < t £ M} U {00}. We write H XtV f° r the hyperbolic arc joining x 
and y. A hyperbolic polygon is a polygon composed of hyperbolic arcs. 

Through the course of this paper, when a vertex of a hyperbolic arc is in 
Q it will always be assumed to be in the form # with a, b £ Z, (a, b) = 1 
and b > 0. If the vertex is 00, we will write it either ^ or ^ (depending on 
if it is the leftmost or rightmost element of a Farey sequence) . 

Let p = \ + and let T be the hyperbolic triangle with vertices p, p 2 
and 00. Then T is a fundamental domain for PSL2(Z) ([Kob93] Prop. III.l). 
Let E e be the edge joining i to 00, E Q be the edge joining p to 00, and Ef 
be the edge joining i to p. Then we call an arc A in the upper half plane 
an even edge (resp. odd edge, resp. f-edge) if A = jE e (resp. A = jE Q , 

resp. A = 7-E/) for some 7 G PSL 2 (Z) (See Figure 1). E e and ^ E e 

together form a hyperbolic arc from to 00, and in general even edges come 
in pairs joining rational numbers |, and y with \a'b — ab'\ = 1 because of 
the following lemma: 

Lemma 2. 1/7 6 PSL2(Z) and a\/b\, 02/62? a'i/b'i, and a' 2 /b' 2 are rational 
numbers in simplest form such that 

7(01/61) = ai/6'i, and 7(a 2 /6 2 ) = a 2 /6 2 



a 2 &i — ai&2 = ct^b'i — a'ib'2 



1 



CHRIS A. KURTH AND LING LONG 



Figure 1. Even edges are thin, odd edges are thick, and 
f-edges are dashed 



Proof. If 7 = then 

,a\. Aai+Bbi a\ .a 2 . Aa 2 + Bb 2 

T(tt) = 7^— -Fa- = 77 and 7(t-) 



6i ' Cai + Dh b[ Iy b 2 ' Ca 2 + Db 2 b' 2 ' 
So: 

ai^-a' 2 6i = (Aa 1 + Bb 1 )(Ca 2 + Db 2 )-(Aa 2 + Bb 2 )(Ca 1 + Db 1 ) 
= ADa\b 2 + BCa 2 bi - ADa 2 b\ - BCa\b 2 
= (AD - BC)( ai b 2 - a 2 h) 
= (aib 2 - a 2 bi) 

□ 

So the quantity a 2 b\ — a\b 2 is invariant under transformations in Y"SL 2 (7L). 
Note that even edges, odd edges and free edges only map to even edges, odd 
edges and free edges respectively under transformations 7 € PSL<2(Z). 

Definition 2. A special polygon P is a convex hyperbolic polygon together 
with a side pairing defined in the following way: The polygon is such that: 

(1) The boundary of P consists of even and odd edges. 

(2) The even edges of P come in pairs, each pair forming a hyperbolic 
arc between elements of Q U {00}. 

(3) The odd edges of P come in pairs, each pair meeting a vertex with 
initial angle 

The sides of the polygon are denoted as follows: 

(1) Each odd edge is called an odd side. 

(2) As even edges come in pairs, either each edge of the pair is an even 
side, or the union of the two edges (a semicircle) is called a free 
side. 

The side pairing on the edges is defined as follows: 
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(1) Each odd side is paired with the odd side it meets at an angle of 
This is called an odd pairing. 

(2) Each even side is paired with the even side with which it forms a 
semicircular arc. This is called an even pairing. 

(3) There are an even number of free sides and they are partitioned into 
sets of two, each called a free pairing. 

We will always assume that and oo are vertices of P. 

The sides of a special polygon P have a natural orientation obtained 
by tracing the perimeter of the polygon in a certain direction. If {s, s'} 
is a side pairing then there is a unique 7 € PSL<2(Z) such that 7 maps s 
to s' in an orientation-reversing manner. We call this the side pairing 
transformation associated with the side pairing, and we let Tp be the 
group generated by all the side pairing transformations of P. Note that it 
doesn't matter which side we pick for s and which for s' because the two 
possible 7's are inverses of each other. Also note that if s is an even side 
(resp. odd side) then 7 is order 2 (resp. order 3). 

Two theorems of Kulkarni are fundamental here: 

Theorem 1 ([Kul91] Theorem 3.2). If P is a special polygon then P is 
a fundamental domain for Tp. Moreover, the side pairing transformations 
{7i} are an independent set of generators ofTp (i.e. the only relations on 
the 7j 's are jf = 1 or 7? = 1 for any finite- order Ji's). 

Theorem 2 ([Kul91] Theorem 3.3). For every T C PSL^Z) of finite index, 
there is a special polygon P such that T = Tp . 

Proof. [Kul91] and also follows from the proof of the algorithm in Section 
4. □ 

Note that although it is true that any subgroup of PSL2(Z) with funda- 
mental domain F is generated by the transformations that map its edges 
together, the fact that the set of generators of a special polygon is an inde- 
pendent set of generators is something special to the special polygon. For 
example T(2) has a fundamental domain shown in Figure 2. There are 

six sides, and the three side pairing transformations are ^ ' f 2 1 

/3 -2\ 

and ( ^ -j) But this is not a independent list of generators because 
3 -2V 1 fl 2\ fl 



2 -\) \0 I J \2 1, 

A special polygon for T(2) is shown in Figure 3. The pairing transfor- 
mations from the special polygon are ^ ^ and ^ ^ • These are 
independent generators of T(2). 

3.2. Farey Symbols. Recall that the classical Farey sequences F n are con- 
structed by taking all the rational numbers < a/b < 1 with denominator 
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Figure 3. A special polygon for T(2) 

at most n and (a, b) = 1 and writing them as a finite sequence in ascending 
order {ao/fro, ■ ■ ■ , a n /b n }. Then for each i we have aj+i&i — aj&i+i = 1. We 
are interested in sequences that satisfy this condition. 

Definition 3. A generalized Farey sequence is a finite sequence: 

III l\ 
\ q > • • • j x m q J 

such that: 

(1) Each Xi = ai/bi is a rational number in reduced form with 6j > 0. 
Additionally, we often consider x_i = ^ and x n+ \ = g. 

(2) // we Zei a_i = —1, b-i = 0, a n +i = 1, and 6 n +i = then 

a i+1 bi - aib i+1 = 1 (1) 

for —l<i<n. 

Note that this definition forces xo and x n to be integers. We will always 
assume X{ = for some i. 
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Definition 4. A Farey symbol is a generalized Farey sequence with some 
additional pairing information. Namely, between each adjacent entries Xi-i 
and Xi we assign a pairing pi which is either a positive integer called a free 
pairing or the symbol "o" called an even pairing or "•" called an odd 
pairing. Each integer that appears as a free pairing appears exactly twice 
in the pairing information. 

So if P is a special polygon, let xq, . . . ,x n be the vertices of P lying in 
Q listed in ascending order. Recall these vertices satisfy aj+i&j — aj&j+i = 
1. Then {-^p, xq, . . . , x n , g} is a generalized Farey sequence. We make a 
Farey symbol out of the generalized Farey sequence by adding the pairing 
information in the obvious way. 

On the other hand, if F is a Farey symbol we can construct a special 
polygon for F. For adjacent entries of the Farey sequence and x«, if pi 
is a free pairing or an even pairing we let P have as a side the hyperbolic arc 
joining X{-\ and X{. Otherwise if it is odd we let 7 be the unique element 
of PSL^Z) such that 7(0) = X{-\ and 7 = X{ and join and Xi by the 
arcs j(Hq jP ) and j(H Pj00 ). Thus we get a hyperbolic polygon which is made 
into a special polygon by adding pairing information in the obvious way. 

Example 1. T(2) has a Farey symbol —00 . j \ | 00 

1221 

3.3. Generators. If P is a special polygon for a group T then V is inde- 
pendently generated by the transformations mapping each side to its paired 
side. If F is a Farey symbol: 

-00 a /b ai/bi • • • a„_i/b n _i a n /b n 00 

PO Pl P'2 Pn-1 Pn Pn + 1 

then we can explicitly give formulas for the 7 corresponding to a given side 
pairing. 

Theorem 3. Suppose (ai/bi, dj+i/ftj+i) are two adjacent vertices of F. 
Then if the pairing between them pi + \ is an even pairing, let: 

'a i+1 b i+1 + aA -af - af +1 
b l + b l+i -a i+ ib i+1 - aibi 
If Pi+i is an odd pairing, let: 



a-i+ih+i + aib i+ i + aibi -af - aja i+ i - af +1 
b\ + hb i+ i + 6f +1 -a i+1 b i+1 - a i+1 bi - ajb, 



And if Pi+i is a free pairing that is paired with the side between a^/bk and 
att+i/bk+i, let: 

q _ ( ak+ih+i + akh —akOi — ak+\a i+ i 
%+1 \bkbi + bk+ibi + i —ai+ibk+i — aibk / 

Then Gj+i is the side transformation corresponding to the pairing pi + \ . 
Proof. [Kul91] Theorem 6.1 □ 
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3.4. Group Invariants. Several invariants of the group T can be read off 
from the Farey symbol F. Firstly, the number of inequivalent order-2 (resp. 
order-3) elliptic points, e 2 (resp. 63), is the number of even (resp. odd) 
pairings in F. Also, the number of free pairings in F (half the number of 
free edges) is equal to r, the rank of 7ri(r\H) (the fundamental group of the 
uncompactified modular curve). 

To discuss the cusps of T, note that if (x«, Xj+i) is an edge with an even or 
odd pairing, then Xi and Xj + i are equivalent cusps (since Gj+i € T maps Xj 
to Xj+i). Likewise, if (xj,Xj + i) and (xj,xj + i) are paired edges then Xj and 
Xj+i are equivalent cusps and Xj and Xj+i are equivalent cusps. This defines 
an equivalence relation on the vertices of P. The equivalence classes are easy 
to compute, because the defining equivalences occur in a cyclic patten. So 
the number of cusps t can be counted as the number of equivalence classes. 

For an edge (f,gi) let 7= J^) ■ So j' 1 ^) = 00 and 7 "W) = 

0. Then define the width of a vertex x« to be the "width" of 7P at 00. That 
is: 

{l^i-i&i+i — cij+ibj_i| if Xi is adjacent to no odd edge 

|cij_i&j + i — aj + i&j_i| + 1/2 if Xj is adjacent to 1 odd edge 
\ai-\bi + \ — aj + i&j_i| + 1 if Xj is adjacent to 2 odd edges 

The cusp width of a cusp x of V is then the sum of the widths of the vertices 
of P T-equivalent to x. 

r\H is a genus g orientable surface with t points missing, one for each 
cusp. The rank of its fundamental group is r = 2g + t — l, so we can calculate 
the genus g = r ~^ +1 . Moreover, using the Hurwitz formula ([Shi71] Prop. 
1.40) we get the index of T in PSL 2 (Z), n = 3e 2 +4e 3 + 12#+6£-12. An even 
simpler formula for the index comes from noting that n + 2 = 2r + e 2 + &z 
where n + 1 is as in Definition 3. This, combined with the previous formula, 
implies /x = 3n + e^. 

4. Coset Permutation Representation of a Group 

Another method of representing groups that will be useful to us in de- 
termining if a group is congruence is the coset permutation representation 
developed by Millington [Mil69a] [Mil69b]. Let T be a subgroup of PSL 2 (Z) 
with [PSL 2 (Z) : T] = /i and PSL 2 (Z) = U^ =1 ajF a coset decomposition with 
oi\ = I. Let F be the standard fundamental domain for PSL 2 (Z). Then 
uf =1 a i _1 F is a fundamental domain for T. Let 

~*=(-°i «)■ v =(-i £ =G !)• fl ^0 !) 

E and V generate PSL 2 (Z), as do L and R. The conversions between 
them are: 

E = LR~ 1 L, V = R- l L (2) 

L = EV' 1 , R = EV- 2 (3) 
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We have E 2 = V 3 = 1. In fact it is well-known that PSL^Z) is isomorphic 
to the group (cf. [Ran77]): 

PSL 2 (Z) 9* (e, v : e 2 = v 3 = 1} (4) 

For each 7 in PSL2(Z), left multiplication acts on the left cosets of T in 
PSL2(Z) by permutation, i.e. there is a homomorphism 4> '■ PSL^Z) — > S n , 
such that if ^(7) = cr 7 then 70^ = a a r^T. In this way every finite-index 
subgroup of PSL2(Z) is associated with a pair of permutations e = <p(E) 
and v = (p(V) with e 2 = i> 3 = 1 which generate a transitive permuta- 
tion group (transitivity comes from E and V generating PSL^Z)). We 
call (e,v) a coset permutation representation of T and (7,r) an LR- 
representation of T, where I = <f(L) and r = <f(R). Each form can be 
obtained from the other form by the equations (2). Note that 7 G PSL^Z) 
is in T if and only if 7r = T, i.e. c 7 (l) = 1. 

On the other hand, suppose e and v are a pair of permutations on p letters 
with e 2 = v 3 = 1 that generate a transitive permutation group S (such a 
permutation we call valid). Define a homomorphism 99 : PSL^Z) — > S 1 such 
that tp(E) = e and v?(U) = v (This is well-defined because of (4)). Let T = 
{7 G PSL2(Z) : 9?(7)(1) = 1}. Then T is an index-fi subgroup of PSL2(Z). 
Thus we have a correlation between valid pairs of permutations and finite- 
index subgroups of PSL2(Z). To test if A G PSL^Z) is in T we write A as a 
word in L and i? (Using, essentially, the Euclidean Algorithm) and replace 
L and R with the permutations / and r. If the resulting permutation fixes 
1 then A is in T. 

If one of the cosets is fixed by e, say e(i) = i, it corresponds to an elliptic 
element in T, for Ea^T = otiT means oT l EaiT = T, meaning a^ 1 Eai (which 
is order 2) is in T. So e2, the number of inequivalent elliptic elements of 
order 2 in T, is equal to the number of elements fixed by e. Similarly, is 
the number of elements fixed by v. 

The cusp width of T at 00 is the smallest positive integer n such that 
L" £ T. Thus the cusp width at infinity is the order of the cycle in (p(L) 
which contains "1". Likewise, suppose i is in a cycle of length k in ip(L), 
i.e. L k ai T = ai F, but L n ai F ^ ai T for < n < k. Then a~ x L k a.i G T, 
but a^ x L n ai ^ T for < n < k. If q = a^ 1 ^ then a~[ l L k aiq = q but 
a^ 1 L n a.iq 7^ g for < n < fc. Thus (x[ l L k oii is a generator for the stabilizer 
of the cusp q, and this cusp has width k. 

5. Algorithms 

5.1. Calculating a Farey Symbol. Recall that T is the standard fun- 
damental domain for PSL2(Z), and let T* be the hyperbolic triangle with 
vertices p, i and 00 (So T = T* U (-T*)). & = {7T : 7 G PSL 2 (Z)} is 
a tessellation of the upper half plane and any finite index subgroup V has 
a fundamental domain which is a simply connected union of ^-tiles. Let 
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Figure 4. A hyperbolic triangle 



ZT* = { 7 T* : 7 € PSL 2 (Z)} U h(-T*) : 7 € PSL 2 (Z)}. ^* is also a tessel- 
lation of the upper half plane, and we will construct a fundamental domain 
for T out of ^*-tiles. The starting point for our construction will be the six 
tiles around an odd vertex. The following lemma shows this is a reasonable 
starting point: 

Lemma 3. Let T be a subgroup of PSL^Z) with index > 3. Then the 
stabilizer of p = \ + ^-i or p — 1 = — ^ + ^-i is trivial (i.e. one of these 
points is not elliptic in T). 

Proof. If the two stabilizers are not trivial then they must be T p = { /, A, A 2 } 

and r p _i = {I,B,B 2 } where A = ^ and B = jV But A 

and B generate an index-2 subgroup of PSL2(Z). So T is either the (unique) 
index-2 subgroup of PSL^Z) or PSL2(Z) itself. And if the index of T in 
PSL/2(Z) is bigger than 2, at least one of A and B cannot be in T. □ 

So if T is not PSL<2(Z) or T2, the unique index 2 subgroup of PSL<2(Z) (cf. 
[Ran77]), then the hyperbolic triangle with vertices either 0, 1 and 00, or 
— 1,0 and 00 (cf. Figure 4) is contained in a fundamental domain of T. The 
triangle is made of 6 j7*-tiles. We will make a polygon P starting with this 
triangle, then attach ^*-tiles to P and assign partial pairing information 
to sides until we get a fundamental domain for T (at which point all the 
pairing information will be filled in). In the algorithm we will say a ^*-tile 
T is adjoinable to P if T is adjacent to a tile of P and if P U T is contained 
in some fundamental domain of V. Note that if T is adjacent to P with 
adjacency edge e and if e cannot be paired with any other edge of P then 
T is adjoinable. 

Algorithm: 

(1) If T = PSL<2(Z) let P be the special polygon with Farey symbol 



—00 



0. 



00 



or if T = T2 let P be the special polygon with Farey symbol 
—00 00 
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Figure 5 



In either case return P and terminate. 

(2) If ^ | ^ is not in T then let P be the hyperbolic polygon with 

vertices 0, 1, and oo. Otherwise let P be the hyperbolic polygon 
with edges —1, and oo. 

(3) If any of the three sides of P map to each other by a 7 £ T, assign 
that pairing to the side. (Note that initially all sides are even sides). 

(4) P is now a polygon where every side is either: 

(a) even and already paired. 

(b) odd and already paired. 

(c) even and unpaired. 

(5) Pick an unpaired even side e. Figure 5 shows the typical case (The 
other cases are the same as this case with everything translated 
by some 7 G PSL^Z)). Since e is unpaired, T\ and T2 must be 
adjoinable. If o\ and 02 are the new odd edges of P after adding T\ 
and T2 to P then either 701 = 02 for some 7 6 T, or there is no such 
7. If there is 7 pair the two edges and go to Step 3. 

(6) If o\ doesn't pair with 02 then it doesn't pair with any other side 
because the only other unpaired sides are odd. So tiles T3 and like- 
wise T4 are adjoinable. Each of these tiles has a free edge and the 
free edges cannot pair with each other (because their common vertex 
would have an internal angle of 4^, so the pairing transformations 
would make things overlap), so T5 and Tg are adjoinable. 

(7) We've now added 6 ^*-tiles to P (One even triangle). If either of 
the new even edges pair with any of the old unpaired even edges 
then assign that pairing. 
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(8) If all the sides of P are paired then we are done. Otherwise go to 
Step 4. 

The output of the algorithm is a special polygon P with Tp = V. Note 
that the algorithm must terminate, because a fundamental domain of V 
has hyperbolic area |[PSL/2(Z) : T] and a single ^*-tile has area |. So 
for P to be contained in a fundamental domain of T it can have at most 
2 • [PSL 2 (Z) : r] Utiles. 

To effectively implement the algorithm we use Farey symbols. We need 
only a way to test for group membership. Note that iipi/qi and pi + \/qi + \ 
are two adjacent vertices of the fundamental polygon then the hyperbolic 
triangle added to the edge H XijXi+1 in Step 6 is the triangle with vertices 
Pi/q-h Pi+i/qi+i, and (pi +p i+ i)/(qi + q i+ i). 

So given a finite-index subgroup T of PSL^Z), if we have a way to test 
for group membership we can calculate a Farey symbol by the following 
algorithm: 

Algorithm for calculating a Farey Symbol: 

(1) If ^ |^ and ^ are in T then T = PSL<2(Z), so return 

— oo , oo 

o • 

and terminate. If ^ ^ ^^j and ^ j are i R ^ then r = T2, 
so return 

—00 ^ 00 

• • 

and terminate. 

(2) If ^ ^ T then let F be the (partial) Farey symbol: 



-00 j j 00 



Otherwise let F be: 

—00 ^ p ^ ^ ^ OO 

(3) For each i with < i < n + 1, if the pairing between and Xj is 
not filled in then check if it can be paired with itself (even or odd 
pairing), or if it can be paired with another unpaired edge (i.e., check 
if the appropriate Gi is in T). Wherever something can be paired, 
assign that pairing. 

(4) If all edges are now paired, return F and terminate. 

(5) If there is still an unpaired edge, say between Pi/qi and 

make a new vertex (pi +Pi+i)/ (qi+qi+i) with no pairing information 
on the edges adjacent to it. Go to Step 3. 
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The output is a Farey symbol for T. 

5.2. Group Membership. The following algorithm described in [LLT95] 
tests if A £ PSL2(Z) is an element of the group corresponding to a Farey 
symbol F. We will need a lemma about even lines: 

Lemma 4. Let I be an even line ( a semicircle on the upper half plane with 
rational endpoints a/b and a' /b' such that \ab' — a'b\ = 1). Let P be a special 
polygon in HI. Then either I C P or If] P = 0. 

Proof. [LLT95] Proposition 2.1 □ 
Let r be a finite index subgroup of PSL2 (Z) and A an element of PSL2 (Z) . 

a -{z i) <5 » 

A maps the even line Ho yOQ to I = H c y d ^ C{j / do . By the lemma, either I C P 
or it is disjoint from P (except possibly at endpoints). If it is disjoint there 
is an edge which it is naturally "closest" to (In a sense discussed in [LLT95] ) . 
The idea of the algorithm is to translate P across the "closest" edge until 
P intersects H c i d c i i d i , at which point A will be in T if and only if I is the 
image of (0, 00) or the an edge paired with (0, 00). In the actual algorithm 
we work in the other direction, translating the even line instead of the special 
polygon. 

Algorithm: [LLT95] 

Let k = and F be a Farey symbol for V with as one of its vertices. 
Without loss of generality, we can assume % < 

(*k a k 

c! 

(1) There are two possibilities: If and ^ are both vertices of P then 

terminate. Otherwise we must have Xi < < ^ < Xj+i with at 

least one "<" a strict inequality. 

(2) Let gi + i be the generator corresponding to the pairing pi + \ (recall 
this is the transformation mapping I = H c > i d i tCk /d k to its paired 
side). If Pi+\ is a free or even pairing, let a& = gi+i- If Pi+i is an 
odd pairing, let m = where x« = 77- , Xj+i = r 1 ^. Then the 

interval (^ L , %) must be between either xi and m or between m and 

k k 

x i+1 . If ^ < m, let a k = g i+1 . Otherwise let a k = g^ v 

(3) Let ^±1 =a k -% L , ^ = a k -^- Replace k with k + 1 and go to 

a.k+1 a k a k+1 a k 

Step 1. 

c' 

The algorithm returns and which are two vertices of P, and a list 
of aj's. 

Theorem 4. TTie algorithm terminates, and A is in T if and only if one of 
the following is true: 
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and and oo are adjacent vertices with an 



even pairing between them. 
Proof. [LLT95] 



□ 



In addition, if A is in T, A can be written as a word in the generators of 



one of the generators for F. 

5.3. Coset Representatives. Let T be a group with special polygon P. 
Let T be the hyperbolic triangle with vertices i, p, and oo. By the construc- 
tion of P, T is contained in P. The set of 7 <G T such that 7T is in P is a 
set of coset representatives of F. 

Let a,i/bi and ai + \/bi + i be a vertex of the special polygon, and let T = 



0. Let Wi be |aj_i6j + i — aj+iA-i| if the pairing between Oijbi and ai + \/bi + \ is 
not an odd pairing and laj-iA+i — aj+i6j_i| + l if it is. Then Wi is the number 
of ^*-tiles of the form 7T in P. Thus a list of left coset representatives for 



5.4. Congruence Testing. Let T be a finite index subgroup of PSL2(Z). 
Lang, Lim and Tan give a test purely in terms of Farey symbols to determine 
if r is a congruence group [LLT95]. Their test relies on Wohlfahrt's Theorem 
[Woh64] which says that if T has level N then T is a congruence group if and 
only if T contains T(N). In Lang, Lim and Tan's test, if T has level N one 
computes a Farey symbol for T(N), giving a complete set of generators for 
r(7V). One then checks if each of these generators is contained in T using 
the above algorithm. The difficulty with this algorithm is that the index of 
T(N) increases very quickly with N, so if T has large level, the calculation of 
a Farey symbol for F(N) can be very lengthy, even if T has relatively small 
index. 

Another test for congruence was developed by Tim Hsu using Millington's 
coset permutation representations [Hsu96]. If we have an LR-representation 
of r there is a list of relations that are satisfied if and only T is congruence. 

To calculate an LR-representation from a Farey symbol, use the above al- 
gorithm to calculate a list of left coset representatives «j € PSL2(Z) where 
PSL2(Z) = Uf=i a «r. To calculate for instance, recall that / is the per- 
mutation such that LctiT = a^T. So / sends i to the unique j such that 
aJ 1 Lai £ T. So we run through every 1 < i < fi and calculate the permu- 
tation, r can be calculated similarly. (Actually, although we need / and r 





Then Lp 1 (a i /b i ) = 00 and 92 1 (a i+ i/b i+1 ) = 



ris ur=o{ T ^V:0<j<^}. 
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it is easier to calculate e and v, because we know beforehand that they are 
order 2 and 3 respectively. Then I = ev' 1 and r = ev~ 2 ). 

Knowing I and r we can directly apply Tim Hsu's congruence algorithm 
[Hsu96]. Depending on the order of I, (i.e. the level of T) there are different 
lists of relations of I and r that are satisfied if and only if T is congruence. 
For example, if N is the order of / and N is odd then T is a congruence 
group if and only if r 2 l~? is the identity permutation (where ^ is the inverse 
of 2 modulo N). 

6. Implementation 

Helena Verrill has a written a MAGMA package for working with Farey 
symbols for congruence groups. Also, for congruence or noncongruence 
groups, the algorithms described above have been implemented by the first 
author as a collection of functions for SAGE. The package and basic exam- 
ples may be downloaded at: 

http : //www. public . iastate . edu/~kurthc/research/ index .html 
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